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ABSTRACT 

Determining the electric field distribution on the Sun's photosphere is es- 
sential for quantitative studies of how energy fiows from the Sun's photosphere, 
through the corona, and into the heliosphere. This electric field also provides 
valuable input for data-driven models of the solar atmosphere and the Sun-Earth 
system. We show how observed vector magnetogram time series can be used to 
estimate the photospheric electric field. Our method uses a "poloidal-toroidal 
decomposition" (PTD) of the time derivative of the vector magnetic field. These 
solutions provide an electric field whose curl obeys all three components of Fara- 
day's Law. The PTD solutions are not unique; the gradient of a scalar potential 
can be added to the PTD electric field without affecting consistency with Fara- 
day's Law. We then present an iterative technique to determine a potential 
function consistent with ideal MHD evolution; but this field is also not a unique 
solution to Faraday's Law. Finally, we explore a variational approach that mini- 
mizes an energy functional to determine a unique electric field, a generalization 
of Longcope's "Minimum Energy Fit" . The PTD technique, the iterative tech- 
nique, and the variational technique are used to estimate electric fields from a 
pair of synthetic vector magnetograms taken from an MHD simulation; and these 
fields are compared with the simulation's known electric fields. The PTD and 
iteration techniques compare favorably to results from existing velocity inversion 
techniques. These three techniques are then applied to a pair of vector magne- 
tograms of solar active region NOAA AR8210, to demonstrate the methods with 
real data. 

Careful examination of the results from all three methods indicates that evo- 
lution of the magnetic vector by itself does not provide enough information to 
determine the true electric field in the photosphere. Either more information 
from other measurements, or physical constraints other than those considered 
here are necessary to find the true electric field. However, we show it is possible 
to construct physically reasonable electric field distributions whose curl matches 
the evolution of all three components of B. We also show that the horizontal 
and vertical Poynting fiux patterns derived from the three techniques are similar 
to one another for the cases investigated. 
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Introduction 



(A version of this manuscript with higher quahty images is at http : //tinyurl . com/yzxn922 ). 



The availabihty of frequent, high quahty p hotospheric vector m agnetogram observations 
from ground-based instrumen ts such as SOLIS (|Hennev et aLll2009[ ). space-based instruments 
such as the Hinode/SOT SP (ITsuneta et aLll2008[) and the planned HMI instrument on SDO 
( Scherrer fc The HMI Team 20051 ) motivates a fresh look at how these data can be used for 
quantitative studies of the dynamic solar magnetic field. 

While the reduction of the observed polarimetry data into maps of the three magnetic 
field components is itself a challenging problem, we assume for simplicity in this paper that 
that problem has been solved, and that time sequences of error-free vector magnetic field 
maps have been obtained over some arbitrary field of view at the solar photosphere. We will 
not address the implications of errors in the measurements, problems with data reduction, 
or uncertainties such as the resolution of the 180 degree ambiguity. 

There are many possible uses of vector magnetic field maps of the photosphere. We 
will focus on just one: the use of time sequences of vector magnetograms to determine the 
surface distribution of the electric field vector on the Sun. 

Attempts to measure electric fields on the Sun have been made using sp ectropolarimetric 
techn iques designed to measure the linear Stark effect in H I spectral lines ( iMoran fc Foukal 
199ll ). Attempts to measure the electric field in prominences showed no signal above the 
measurement threshold, but these te chniques applied to a small flare surge did result in 
a measurement above the threshold ( iFoukal &: Behrl Il995l ). In this paper, we attempt to 
determine the electric field in the solar atmosphere indirectly from vector magnetic field 
measurements, using Faraday's Law, rather than appealing to the Stark effect. 

There are several reasons why determining the electric field on the solar surface is useful. 
The simultaneous determination of three-component magnetic and electric field vectors will 
allow us to estimate the Poynting flux of electromagnetic energy entering the corona, as 
well as the flux of relative magnetic helicity. In addition, under the assumptions of ideal 
MHD, where E = — v/c x B, it allows us to estimate the three-component flow field in the 
photosphere. Knowing either the flows or an electric field consistent with Faraday's Law will 
enable the driving of MHD models of the solar atmosphere that are consistent with observed 
data. This is a key requirement for predictive, physics-based models of the solar atmosphere 
that might be used in forecasting applications. 

Most recent research on deriving electric fields in the solar atmosphere has been done 
by either explicitly or implicitly invoking the ideal MHD assumption, E = —v/c x B, 
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Fig. 1. — The left panel shows a vector magnetogram of Active Regio n 8210 taken with th e 
University of Hawaii's IVM instrument. This dataset is described in IWelsch et al.l ( I2004J ). 
Arrows show the directions and amplitudes of and By, and the background image shows 
the amplitude of B^. The right panel shows a "vector electrogram" (a three-dimensional 
vector electric field map) of the active region, using the time evolution of B to estimate E. 
Arrows show estimated directions and amplitudes of and Ey, while the background image 
shows the estimated amplitude of Ez- The example shown here displays E computed using 
the "variational" technique ( §3.3p . A detailed discussion of the calculation is in §11 



-4- 



and has focused on deriving two- or three-component flow fields from time sequences of 
magnetograms. Such techniques can be divided into two classes, which we call "tracking 
methods", and "inductive methods". 

Tra cking methods, such as the local correlation tracking (LCT) approach, first devel- 
oped by [November fc SimonI (119881 ). find a velocity vector by computing a cross-correlation 
function that depends on the shift between sub-images or tiles when comparing two im- 
ages. The shift that maximizes the cross-correlation function (or alternatively, minimizes 
an error function) is then taken to be the displacement between the two sub-images; this 
displacement divided by the time between the images is then defined as the average local 
velocity. The velocity field over the entire image is built up by repeating this process for 
all image locations. The LCT technique suffers from two shortcomings: (1) the technique 
does not assume any physical conservation laws, meaning that the derived velocity field may 
not have a physical connection to the real flow field in the solar atmosphere; and (2) the 
technique is intrinsically two-dimensional, and does not account for vertical flows or evolving 
three-dimensional structures in the solar atmosphere. On the other hand, LCT techniques 
offer the advantage of being simple and robust, and are able to use non-magnetic data, such 
as white-light or G-band images for estimating flow fields, though the results must then be 
interpreted carefully. Examples of LC T techniques in curre i it use in the so lar research com- 



muni ty include the implementation bv [November fc SimonI ( 19881). FLCT (Fisher &: Welsch 



20081), and Lockheed-Martin's LCT code (ITitle et al. 



1995 



Hurlburt et al 



19951 ) 



Inductive methods of flow inversion from magnetograms were pioneered by lKusano et al. 



(120021 ). who used a combination of horizontal flow velocities derived from LCT techniques 
applied to the normal component of vector magnetograms, along with a solution to the 
vertical component of the magnetic induction equation, to derive a three-component flow 
field from a sequence of vector r nagnetograms. An a l ternat e approach, using the same idea, 
but adding the interpretation of Demoulin fc Berger (l2003h plus a Helni ho ltz decomposition 
of the "flux transport velocity" was proposed by I Welsch et al.l (120041 ). iLongcopd (120041 ) 
combined the vertical component of the induction equation with a variational constraint 
that minimizes the total kinetic energy of the photosphere while still obeying the normal 
component of the induction equation. Additional techniques minimize a localized error 
functional wh i le ens uring that the vertical component of the induction equation is satisfied 
JSchucklbood 120081 ). 



Kusano et al.l (120021 ) observed that the dB^/dt and dBy/dt components of the induction 



equation involve vertical derivatives of horizontal electric field terms {dEy/dz and dE^/dz, 
respectively). They noted that these electric field components are unconstrained by single- 
height vector magnetogram sequences. Consequently, while the two horizontal components of 
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the induction equation provide more information about the evolution of the vector magnetic 
field, they introduce two more unknowns to the system. 

The primary goal of this paper is to explore the extent to which an electric field consistent 
with the evolution of all three components of B can be derived from a sequence of single-height 
vector magnetograms, despite incomplete information about dE^/dz and dEy/dz. We first 
show that it is possible to derive an electric field whose curl is equal to the time derivatives 
of all three components of B. However, this electric field is not unique, and knowledge of 
additional physics of the electric field formation is necessary to further constrain the solution. 
We then show that a solution for the electric field, under the assumptions of ideal MHD, 
can be derived through the use of an iterative procedure, or alternatively, as the solution of 
a variational problem. We will compare and contrast these solutions with results from an 
MHD simulation, where the electric field is known. To illustrate this methodology with real 
data, we then apply these techniques to a set of vector magnetograms of NOAA AR8210 
from 1 May 1998. Figure [1] shows a vector electric field map (a "vector electrogram") of 
NOAA AR8210, along with a vector magnetogram of the same active region, derived using 
one of the techniques discussed in this paper. 

The remainder of this paper is structured as follows. In ^ we derive solutions to the 
electric field given the vector mag netic field evolution i n a single ho rizontal plane. These 



solutions use the poloidal-toroidal fl Chandr asekharl 1 1 96 ll : lMoffattlll978l ) decomposition of the 
electric field. 

In ^ we show how additional physics describing the electric field can be included by 
exploring two approaches, one iterative and the other variational, both aiming for the con- 
struction of an ideal MHD electric field consistent with E = — v/c x B. While the ideal 
MHD assumption is believed to be a good approximation in the photosphere, the variational 
formalism has been generalized to include non-ideal effects. In §3] we also compare and con- 
strast electric fields derived with our techniques to a test case where the true solutions are 
known, and in §1] with an example using real vector magnetic field data. We discuss and 
summarize our results in 35] 



2. Poloidal- Toroidal Decomposition 
2.1. Decomposing the Magnetic Field and its Time Derivative 



The poloidal-toroidal decomposition (henceforth PT D) of the mag netic field into two 
scalar potentials is well-known among dynamo theorists (jMoffattlll978l ) and has been used 
extensively in MHD models of the solar interior that employ the anelastic approximation 
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appears to have been introduced by IChandrasekharl (jl96ll ). Here, we will briefly discuss the 



PTD of the magnetic fleld, but will focus most of our effort on the PTD of the partial time 
derivative of the magnetic fleld, because of its connection to the electric fleld. Here we use 
B to refer to a snapshot of the magnetic fleld within the vector magnetogram fleld of view, 
and B to refer to its partial time derivative. 

Given a snapshot of the three-component magnetic fleld distribution in Cartesian coor- 
dinates, one can write B as follows: 



B = V X V X i3z + V X J^z. 



Here, J is referred to as the "toroidal" potential, and B as the "poloidal" potential. The 
vector potential A is then given by 



A = V X i3z + + V^, 

where ^ is a gauge potential, left unspecifled. 



(2) 



Taking the partial time derivative of equation ([T]) leads to an equation of exactly the 
same form for B in terms of the partial time derivatives B and J: 



B = VxVxi3z + VxJ^z 



(3) 



The vector z is assumed to point in the vertical direction, i.e. normal to the photosphere. 
A subscript z will denote vector components in the z direction, and a subscript h will denote 
vector components or derivatives in the locally horizontal directions, parallel to the tangent 
plane of the photosphere. 



It is useful to rewrite equations ([T]) and ([3]) in terms of horizontal and vertical derivatives 



as 



B = V, {^-^ + V;, X Jz - V^Sz, 



and 



dB 

dz 



B = V;, I — 1 + V/, X jz - vim. 



(4) 



(5) 



The PTD of B has a useful connection to observation. Examining only the z-component 
of equation ([5]), one flnds 

-B. , (6) 
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where Bz is the partial time derivative of the vertical component of the magnetic field, and 
is the horizontal contribution to the Laplacian. Thus knowledge of Bz in a layer yields 
a solution for B by solving a horizontal, two-dimensional Poisson equation. 

Taking the curl of equation ([5]) and examining only the z-component of the result, one 

finds 

VlJ = -{47T/c)jz = -z ■ (V X B,.). (7) 

Knowing the time derivative of the horizontal field B/^ in a layer, and hence the vertical 
component of its curl, determines J7 in that layer, from solutions to another Poisson equation. 

Finally, taking the horizontal divergence of equation ([5]) results in 

VlidB/dz) = Vh ■ B,. (8) 

Here, knowing the horizontal divergence of B^ allows one to determine dB/dz, once again 
by solving a two-dimensional Poisson equation. 

It is worth noting an additional implication of equation (|H]). From the solenoidal nature 
of B it follows that 

V, ■ B, = ~dBz/dz . (9) 

Thus equation ^ can be regarded as the partial z-derivative of equation ([6]), yet no depth 
derivatives of the data were needed to evaluate it. We return to this point later. 

To find the PTD of the magnetic field itself rather than its time derivative, the quanti- 
ties B, J , and dB/dz obey exactly the same Poisson equations ([MH]) above, but with all of 
the overdots in the equations removed. This description of B allows for an alternate com- 
putation of potential magnetic fields valid at the vector magnetogram surface (Appendix 
|A|) . in which the divergence of B/i taken from the vector magnetogram can be incorporated 
into the solution. This formulation for the potential magnetic fiel d has applicat ions for: («) 



estimating the flux of free magnetic energy, e.g.^ equation (23) of IWelschI (120061 ). where one 
needs to subtract the measured and the potential-field values of the horizontal magnetic 
field, and (?«) computing the vector-potential Ap of the potential field, useful in estimat- 
ing the magnetic helicity fiux. The horizontal magnetic field can be decomposed into the 
potential-field contribution ^ hidB/dz) and the non-potential contribution B^ = x J'z 
(see Appendix |A]) . Welsch's fiux of free energy then becomes 

S'r = ^cE, X Bf , (10) 
where methods of estimating E will be discussed below. 
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2.2. Finding an Electric Field from Faraday's Law 

Now, compare equations (|3]) and ([5]) with Faraday's law relating the time derivative of 
B to the curl of the electric field: 

B = -cV X E (11) 
Equating the expressions for B, one finds this expression for cV x E: 

cV X E = - V X V X i3z - V X jz (12) 
= -Vh{dB/dz)-Vh^ ji' + ^lm. (13) 

Uncurling equation f|T2|) yields this expression for the electric field E: 

cE = - V X i3z - jTz - cVi) = cE^ - cVi/j. (14) 

In the process of uncurling equation f|T2|) . it is necessary to add the (three-dimensional) 
gradient of an unspecified scalar potential ip to the expression for the total electric field. The 
potential ip can be equated to — ^/c + A, where ^ is the gauge potential from equation 
and A is some other potential function. However, we simply let the electric field potential 
ip absorb both contributions, and will not use the gauge ^. The part of E without the 
contribution from —'Vip will be henceforth denoted E"^, the purely inductive contribution to 
the electric field. 

We have derived an expression for V x E, including the two horizontal components 
of the induction equation, simply by using time derivative information contained within a 
single layer, even though those two components of the induction equation include vertical 
derivatives of and Ey. This was made possible by equation (jS]), which enabled the 
evaluation of the needed depth derivatives through the relation V ■ B = 0. 



2.3. Boundary Conditions 

To solve the three Poisson equations (E]), ([7]), and ([8]), one must consider their boundary 
conditions. For many of the MHD simulation test cases we have used, the simulation vector 
fields obey periodic boundary conditions, which makes the problem straightforward: one 
can use Fast Fourier Transform (FFT) techniques to solve the Poisson equations without 
special consideration for boundary conditions (see Appendix [B|) . For an arbitrary vector 
magnetogram taken over a finite area, however, the magnetic field will generally not be 
periodic, but will be determined by the measured fields on the boundary. For equations 
d?]) and ([H]), the horizontal components of B determine the boundary conditions for J' and 
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dB/dz from the x and y components of the primitive equation ([5]): 

^ _d_dB_^dj_ 

dx dz dy ' 

and 

B = (16) 

^ dy dz dx 

From these equations, coupled Neumann boundary conditions can be derived: 

= (17) 
dn dz ^ ds '' 



and 



^ = -5. + -—, (1^ 

dn ^ ds dz ' 



where subscript n denotes components or derivatives in the direction of the outward normal 
from the boundary of the magnetogram, and subscript s denotes components or derivatives 
in the counter-clockwise direction along the magnetogram boundary. 

The choice of boundary conditions for B in equation ([6]) can have subtle effects on 
the solution Bz- If homogeneous Neumann boundary conditions (derivative normal to the 
magnetogram boundary specified with zero slope) are applied to B, then the horizontal curl 
of around the magnetogram boundary necessarily vanishes, meaning the average value 
of Bz within the magnetogram is forced to zero. If the average value of Bz is nonzero, the 
solutions will not reflect this. This problem can be corrected post facto, however, by adding 
a correction term to the electric field - see Appendix O 

Applying homogenous Dirichlet boundary conditions (setting B to zero at the bound- 
aries) for the solution of equation ([6]) for B can also result in artifacts if care is not taken. 
Homogenous Dirichlet boundary conditions on B imply no net change in B across the mag- 
netogram, resulting in average x— and y— components of the electric field that are zero. 
Some evolutionary patterns, such as the emergence and separation of a simple magnetic 
bipole oriented in the x— direction should result in a non-zero average value of Ey, because 
the opposite polarities have opposite velocities. Thus, if Dirichlet boundary conditions are 
used in the solution of equation for B, some other technique must be used to find the 
magnetogram- averaged values of and Ey. 

The coupled boundary conditions fll7H18l) for the two Poisson equations ([3IH]) are degen- 
erate, in that there is a family of coupled non-zero solutions to the homogeneous Cauchy- 
Riemann equations for J' and dB/dz: 
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and 

dj 




3. 

which satisfy the boundary conditions for zero time derivative of the horizontal magnetic 
field on the boundary. These solutions can be added to solutions for j and dB/dz without 
changing the time derivative of the horizontal field on the boundary. Solutions to equations 
(119112 op also are harmonic, i.e. each of the two solutions also obey the two-dimensional 
(horizontal) Laplace's equation. 

The boundary conditions needed to solve the Poisson equations for i3, J , and dB/dz (if 
one needs to perform the PTD on B, rather than B) are identical to the boundary conditions 
described above, but with all of the overdots removed. The Cauchy Riemann degeneracy 
described above also applies to the solutions of J and dB/dz. 

Even without considering V?/', the degeneracy of the solutions for j and dB/dz for non- 
periodic boundary conditions means that particular solutions of the homogeneous Cauchy- 
Riemann equations can be added to the solutions for the inductive electric field cE^ without 
affecting V x E. This means there is some freedom to specify the solution of J or dB/dz at 
the boundary. This freedom cannot be applied to both functions simultaneously, however, 
since the two solutions are coupled together. In practice, this means that one could choose to 
set the derivatives parallel to the boundary of one of these two functions to zero, for example, 
while still obeying the coupled boundary conditions (fT7|) and ( IT8|) . For an illustration of this 
degeneracy applied to the PTD solutions for B, see Figure [2j 

When periodic boundary conditions cannot be assumed, equation ([H]) shows that the 
electric field on the boundary of the magnetogram determines the boundary conditions for 
B and ip: 

dB d^ 



and 



M = -eE (22) 



These coupled Neumann boundary conditions have a similar form as those for J' and dB/dz, 
(equations [TTlfTS] ) except that they depend on horizontal electric fields rather than on time 
derivatives of the horizontal magnetic fields. The practical difficulty with using equations 
(12111221) is that generally one does not know the behavior of the horizontal electric field vector 
on the boundaries if the boundaries are in regions of strong, evolving magnetic field . On 
the other hand, if one can take the boundaries along regions with weak average magnetic 
field strength, one could probably set En and Es to zero along the boundaries. 
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Fig. 2. — The top two panels show Poisson equation solutions for dB/dz (left) and J (right) 
within the AR 8210 magnetogram field-of-view, under the assumption that dj/ds = along 
the magnetogram boundary. The AR8210 vector magnetogram data is discussed further in 
^ The bottom two panels show Poisson equation solutions for the same two functions 
under the assumption that d/ds{dB/dz) = 0. The difference between the upper and lower 
sets of solutions obey the homogeneous Cauchy-Riemann equations. The two different sets 
of functions yield identical values for the horizontal components of the magnetic field. 
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3. Determining the Scalar Potential 

3.1. The Implications of Not Specifying ip 

We can test our approach by applying equations (fT^ and (JH]) to the magnetic evolution 
sampled from a thin slab of an MHD simulation, in which E and V x E are both known. 
How well do the derived results compare with the known electric field, and with the curl of 
that electric field? 

To address this q uestion, we use the results from the ANMHD simulations described 
in Welsch et al. ( 2007D. which have b een used for several studies of velocity field inversions 



(jWelsch fc Fisher! I2OO8I : ISchuckll2008l ). This simulation models a magnetic bipole emerging 
through a strongly convecting layer. To compute the reconstructed distribution of V x E 
and E using PTD, we use time differences from two consecutive output steps from the MHD 
simulation to find estimates of the time derivative of the three components of B. From these 
three components, the three Poisson equations (I6]|8]) can be solved, subject to boundary 
conditions (I17til8p for dB/dz and J. To solve the three Poisson equations along with these 
coupled boundary conditions, we first use a simple successive over-relaxation method to solve 
the equation for i3, assuming homogenous Neumann (zero gradient) boundary conditions. 
We choose Neumann over Dirichlet boundary conditions because Dirichlet boundary condi- 
tions on B lead to zero average horizontal electric field within the magnetogram (see earlier 
discussion in §2.3p . To solve the two Poisson equations ([7}l8]) along with their coupled Neu- 
mann boundary conditi ons (equation s [T7] and [IB]), we have adapted the Newton-Krylov 



solver from RADMHD ( lAbbettl 120071 ) to solve this elliptic system simultaneously. 



The top panels of Figure [3] show the curl of the electric field obtained from time differ- 
ences in the magnetic field between adjacent snapshots of the magnetic field evolution. The 
middle panels of the Figure show the reconstructed distribution of V x E as obtained from 
equation ( IT^ . The bottom panels show scatter plots between the known and reconstructed 
values of V X E. One can see that the reconstructed components of V x E show good 
agreement with the original values from the MHD code. 

If one ignores the contribution from V?/' in equation f[T^ . which is equivalent to assuming 
that E is equal to E-'^, it is straightforward to compare the reconstructed electric field E-'^ with 
that used in the ANMHD code. The top panels of Figure H] show the actual components of 
the electric field used in the MHD simulation, while the middle panels show the electric field 
components from equation f[T^ approximating E by E^. The bottom panels show scatter 
plots of the reconstructed electric field components as a function of the known values. 



The velocity inversion methods tested in 



Welsch et al, 



( I2OO7I ) also included comparisons 
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between the ANMHD and reconstructed elec tric field component s computed by assuming 



cE = —V X B (see their Figures 11 and 12). IWelsch et al.l (120071 ) showed that the simula- 



tion data for B^, computed by differencing two images in time, matched — cz ■ (V x E), 
computed using simple second-order finite difference formulae to evaluate the spatial deriva- 
tives apphed to the average of the electric fields at the two adjacent times. Hence, it is 
possible to make a direct comparison between the PTD derived electric field values reported 
here and those fro m several velocity in version methods, including LCT techniques. The 



results reported in IWelsch et al.l (120071 ) include only "strong field" pixel locations where 
\Bz\ > 370G, 5% of the maximum field strength in this simulation of subsurface magnetic 
evolution. Accordingly, to make a direct comparison, we compute rank-order correlation 
coefficients between the PTD-derived electric field components with those from ANMHD, 
including only the strong field locations. The rank-order correlation coefficients for the x— , 
y—, and z— components of E are 0.76, 0.87, and 0.94, respectively. This metric compare s 



favorably with all of the techniques shown in Figures 11 and 12 of IWelsch et al.l ( l2007l ). 
showing better correlation coefficients for and Ey than all the methods except MEF, and 
comparable correlation coefficients to MEF. All the methods, including PTD, show good 
correlation coefficients for E^. 

Despite the good performance of PTD compared to the velocity inversion methods 
for strong-field locations, a scatterplot comparison between the original and reconstructed 
electric fields is poor when all points are considered (the bottom panels of Figured]). This 
contrasts with the scatterplot comparison between actual and reconstructed values of V x E 
which showed good agreement. 

Why is this? The problem occurs because of the under-determination of the electric 
field. The PTD formalism guarantees that the electric field will obey Faraday's law, but no 
other information about additional physical mechanisms that would determine the electric 
field has been incorporated. Any additional electric field contribution that can be represented 
by the gradient of a scalar function is left completely unspecified. Evidently setting V^/^ = 
in the PTD solution (equation [T^ is inconsistent with the electric field as specified in the 
ANMHD simulation. Therefore, finding an equation that better describes ^ip is essential 
for providing a better performance of the PTD method. 

Since E ~ — v/c x B in these nearly ideal ANMHD data, E is perpendicular to B. Con- 
sidering only the cE^ = — V x — J'z contributions from the poloidal and toroidal terms in 
equation ([H]), and decomposing the vectors into the directions parallel and perpendicular to 
B, one finds contributions that are roughly equal in the parallel and perpendicular directions. 
To better reconstruct the actual electric field, it will be necessary to add a potential electric 
field that largely cancels out the components of E that are parallel to B. The challenge is 




Fig. 3. — Top panels show the values of the x—, y—, and z— components of cV x E = 
—d'Q/dt taken from the ANMHD simulations described in the text. The middle panels show 
the same components of cV x E determined with the PTD formalism (equations [T2lfT3] ). The 
bottom panels show scatter-plots of the original versus derived values of the x— , y— , and z— 
components of cV x E. Note the scatter-plots for the x— and y— components of cV x E are 
not as tight as the z— component. This stems from using finite difference operations for the 
Laplacian that aren't strictly compatible with the finite difference operations we have used 
to evaluate the curl and divergence operators. We have also performed similar comparisons 
for FFT-derived solutions (Appendix |B]), which can be applied here, since ANMHD does 
assume periodic boundary conditions. In that case, the corresponding x,y scatter-plots are 
considerably tighter than the ones shown here. All of the grayscale images are plotted using 
the same scale. 
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Fig. 4. — Top panels show the values of the x— , ,and z— components of cE taken from the 
ANMHD simulations described in the text. The middle panels show the same components of 
cE determined with the PTD formalism (equation [12]) assuming VV' = 0. Bottom panels 
show scatter plots of the x—,y—,z— components of cE derived from the PTD formalism 
versus the ANMHD components of cE. Note that compared to scatter-plots comparing 
V X E between the simulation data and the PTD results, the results here show considerable 
scatter. All of the grayscale images are plotted using the same scale. 
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to derive an equation for from physical or mathematical principles that does this, while 
also yielding a physically reasonable solution for the resulting total electric field. 



3.2. Deriving an Electric Potential I. - An Iterative Approach 

Here, we describe a technique to determine a potential function that is consistent with 
ideal MHD (E ■ B = 0), using a purely ad- hoc iterative approach. The total electric field is 

E = E^ -Vtfj , (23) 

where as before, E''^ is the inductive contribution found using the PTD formalism and —'Vip 
is the potential contribution. We wish to define ip in such a way that the components of E 
parallel to B are minimized. 

In step 1 of the procedure, we decompose V?/' into three orthogonal directions by writing 
Vip = si{x, y)h + S2{x, y)z x b + S3(x, b x (z x b) (24) 

Here, b is the unit vector pointing in the direction of B. 

In step 2, we set si{x,y) equal to E-^ ■ b. This ensures that —'Vip acts to cancel the 
component of E/ parallel to b. The function si{x, y) will remain invariant during the rest of 
the iteration procedure. 

In step 3, given the latest guess for ip, we evaluate the functions S2{x, y) and s^{x, y) by 
dotting each of the vectors on the right hand side of equation f l2^ with equation fl24|) . itself, 
yielding 

z ■ (bfe X Vfe^) 

S2[x,y) = -2 , (25) 

"h 

and 

.3(x,y) = ^ ^ • (26) 

Here, hz and hh represent the amplitudes of b in the vertical and horizontal directions, 
respectively, and hh represents only the horizontal components of b. For the initial guess 
during step 3, the functions S2{x,y) and S3{x,y) are set to zero. 

In step 4, the horizontal divergence of V/j?/' is taken, using the current guesses for 
S2{x, y) and s^lx, y). From equation flM|) this results in the following two-dimensional Poisson 
equation for ip: 

Vl^p = Vh ■ {si{x, y)hh + S2{x, y){z x b) - ss{x, y)Khh) (27) 
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The values for ijj are updated by solving this Poisson equation. 

In step 5, the vertical gradient of ip is updated by evaluating the z-component of equa- 
tion employing the last guess for ss{x,y) (recall that Si{x,y) does not change between 
iterations) : 

si{x,y)b^ + S3{x,y)bf^ . (28) 



dz 

Step 6 consists of evaluating an error term 

|(E/-V^)-b| 
l(E/-VV^)| ^ ' 

If e is sufficiently small, then the iteration sequence can be terminated; otherwise steps 3-6 
are repeated until the sequence converges to the desired error criterion. 

Results of this iteration sequence applied to the PTD solutions, and compared to the 
true electric field results from the ANMHD simulation are shown in Figure O To solve the 
Poisson equation fl7r|) in this example, we take advantage of the known periodic nature of 
the ANMHD solutions, and use FFT techniques to solve the equation for ip. 

From Figure [5l one can see that some of the artificial features in the PTD recovered 
solutions have been improved by applying the iteration scheme described here, such as the 
false bright halo seen in the PTD solution for on the upper left side. Further, scatterplots 
of and Ey comparing results between simulation values and those with PTD only, and 
those with PTD plus the iteration scheme, show clear improvement by adding in the contri- 
butions from this potential function. However, the technique is by no means perfect. After 
applying the correction derived with the iteration technique, other "shadow" artifacts, 
seen as faint vertical stripes below the emerging magnetic field, are visible in e.g. the derived 
map of E^. Computing the rank-order correlation coefficients between the ANMHD and 
iteration method electric fields in the strong magnetic field regions, as we did for the PTD 
solutions, results in values of 0.75, 0.82, and 0.95 for the x—, y— and z— components of E. 
These values are not significantly different than those for the PTD case described earlier. 

In summary, this iteration scheme, or perhaps similar schemes based on related ideas, 
may provide a useful approach for a deriving potential electric field contribution which, when 
added to the PTD solutions, is consistent with ideal MHD. But we must also caution that 
the solutions derived via this method are not unique since the condition E ■ B = does not 
fully constrain the potential function ip. Once a solution has been obtained via this method, 
we can add on any additional solution ip' which obeys the constraint ^ip' -6 = 0, without 
affecting the induction equation or the condition E ■ B = 0. 
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Fig. 5. — Top panels show the values of the x—, y—, and z— components of cE taken from the 
ANMHD simulations described in the text. The middle panels show the same components 
of cE determined with the PTD formalism (equation [12]) adding —'Vijj as determined with 
the "iterative" technique of §3.2[ Bottom panels show scatter plots of the x—, y— ,and z— 
components of cE derived from the PTD plus potential field solutions versus the ANMHD 
components of cE. Note that compared to scatter-plots comparing cE from PTD solutions 
without the potential contribution (Figure [1]), the results here show less scatter. All of the 
grayscale images are plotted using the same scale. 
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It is also not clear whether this iteration technique will converge for all cases, or whether 
the derived solutions are mathematically well posed apart from the uniqueness issue already 
raised. Clearly this area needs further study. 



3.3. Deriving an Electric Potential II. - A Variational Approach 

The dynamics of the solar plasma is determined by the largest forces, which in regions of 
strong magnetic fields will involve Lorentz forces, acting in conjunction with gravity, pressure 
gradients, and inertial terms. To the extent that the electric field is dominated by the ideal 
— v/c X B term, then it is necessary to know the forces acting to determine the velocity field 
in the photosphere to determine the full electric field and thus completely specify ip. 

In §3.11 we demonstrated that vector magnetograms alone contain only partial infor- 
mation about the plasma dynamics - there simply isn't enough information in the magnetic 
field data alone to uniquely specify v or E. Additional information must be obtained either 
from other measurements or by using some other constraint. 

One approach for deriving a constraint equation for is to use a variational princi- 
ple. For example, one could adjust ip such that the electromagnetic field energy density, 
{E^ + 52)/(87r), integrated over the magnetogram, is minimized. Since B itself has already 
been determined by the measurements, this is tantamount to finding ip such that the area 
integral of E"^ is minimized. The motivation for this approach is to reduce or eliminate 
the unphysical electric field "halos" seen in regions of negligible magnetic fields strength in 
Figure m An alternative equation for ip can be derived by minimizing |cE x B/i?^p inte- 
grated over the magnetogram, which is equivalent to minimizing the kinetic energy of fiows 
in the photosph ere, if one assumes E = —v/c x B. This is essentially the approach used by 
Longcope in his MEF technique for deriving flows from magnetograms using only the 



vertical component of the induction equation. 

Here we derive an equation for ip which is sufficiently general that both of these cases 
can be included using the same formalism. We minimize the functional 

L = j dxdyW\x,y)[{Ei - di)/dxf + {El - d^j/dyf + {El - d^j/dzf] , (30) 

where W'^{x,y) is an arbitrary weighting function, and E^, E^, and E!, are the three com- 
ponents of E^ as it defined in equation 



In the photosphere, we believe the electric field E^ — 'Vip is dominated by the ideal term 
-v/c X B, but with a possible additional contribution R, which can represent resistive or 
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any other non-ideal electric field terms. We assume that R is either a known function, is 
determined by observation, or is specified by the user as an Ansatz. Then 

E = -- X B + R = - Vt/- . (31) 

c 

By dotting E^ — with B, one finds 

(E^ - V^) ■ B = R ■ B. (32) 

This equation provides an additional constraint on the potential ijj, allowing us to eliminate 
dip/dz in favor of V/^'?/': 

= B ■ E^ - B;, ■ V;,^ - R- B. (33) 



Note that the functional minimized in equation ( 130|) depends on ip through its depen- 
dence on the total electric field E^ — V^/^. In particular, the z-component of E that appears 
in equation fl30|) depends on — dip/dz. But from the above constraint equation fl33|) we 
can see that 

B,{Ei - ^) = R . B - B, • (E^ - V,V^), (34) 

showing that El = — J cancels out of the variational equation. This result shows that a 
solution of the variational problem for E'^ — V?/^ will be independent of the solutions for 
j and dB/dz as determined from equations ([318]) and boundary conditions (11711181) . Thus 
the fact that dB/dz and j do not have unique solutions does not affect the uniqueness of 
the solutions for the electric field itself: In this approach, changes in J are compensated by 
changes in dip/dz such that Ez is unchanged. 

Performing the Euler-Lagrange minimization of equation (130]) results in a second-order, 
two-dimensional elliptic partial differential equation for ip, 



= V. . ^ V) ( (e; - VWO + ^, \^:'-^-^')--^-^ ) ^ ,35) 



2 



Equation f p5|) involves a combination of both the inductive electric field from equation ([T] 
and the potential contribution from — V/iV) assumin g that E^ is given . Solving this equation 



for ip given E^ is essentially the approach taken by iLongcopd (120041 ) in the development of 
MEF. 



Alternatively, equation ( 135|) can be viewed as a single equation for the sum of both the 
inductive and potential contributions, to be determined simultaneously. If the equation for 
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the total field can be solved, then the potential term can be found afterward by subtracting 
the PTD solution (equation [14]) from the solution for the total electric field. 

Writing the total electric field E as - Vip, or Ez = El- dip jdz, and E/^ = E^ - V/i?/;, 
and noting that equation (IMl) relates to E/^, equation (135!) can be re- written as 

-V;, • (W^V^.) (E X B) X z = 
-z ■ V;, X {W'^IB,) (E X B)^ = 0. (36) 

The variational approach thus leads to a local condition on the quantity (iy^/i?2)(E x B);i, 
namely that it is curl-free, and thus can be represented as the gradient of a two-dimensional 
scalar function. Therefore, we write 

(iyV5.)(cExB), = -V,x- (37) 



We wish to derive an equation for x, starting from equation (1371) and involving only 
known quantities from the vector magnetic field or its time derivatives. The details of the 
derivation are shown in Appendix [Dl The result is 

Qjy /cR B \ / 1 \ / \ 

' + V^- — — z X B J + V^- TT7^(B;, ■ Vnx)^^ = - V^- -^VhX ■ (38) 



dt V B'^ J yW^B"^^ ' J yw^B 

Eventually, we will assume R = 0, but for now we retain it in our formalism so that non-ideal 
effects can be included. 

If one either knows R or sets R ■ B to 0, equation fl55]) is a two-dimensional linear 
elliptic partial differential equation for x, with coefficients that depend on the magnetic field 
components or its time derivatives. Thus we can regard the solution for x as well-defined. 

In Appendix [D] we determined how to find E^ from x via equation (]D3p . To find E/^, 
one can take the cross-product of z with equation (lD2p and use equation (]D3p to find 

cE/. = c^B, - z X V.x - (B. • Vnx) z x B,. (39) 

The Poynting flux S = (c/(47r)) E x B has components that point along the gradient of 
X in the horizontal direction as the definition (!37|) shows. However, we can determine the 
Poynting flux in the vertical direction as well by noting that S ■ B = 0, or B^Sz = —^h ■ ^h- 
Thus we find 
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and 



and 

Assuming that the non-ideal electric field term R is zero or negligible compared to the 
ideal contribution, equation ( l38l) is greatly simplified for the two special cases of = 1/ B'^ 
(minimum kinetic energy) and W"^ = 1 (minimum electric field energy). In those cases, 
equation ( |38|) becomes 

{W^B^ = 1 :) + V, ■ ((B;, ■ VhX)^u) = -^h ■ [bIV ^,x) ■ (42) 

iW^ = 1 :) + V;, ■ ((b, ■ VhX)^h) = -V, ■ [hlVhx] ■ (43) 

In equation fH3|) and represent, respectively, the vertical and horizontal components 
of the unit vector h pointing in the direction of the magnetic field. 

Once the variational equation for x has been solved, how does one relate the total electric 
field to the PTD solutions and the potential field contribution? Since the total electric field 
E = E-^ — V'?/', we can subtract the E^ contribution of equation f|T^ from equations (l39l) 
and ( 1D3I) to derive these expressions for the electric field from V?/': 

R, B B^ 1 

- cV,^ = c-^B, - zxVhX- 1^ {Bh ■ VhX) z X B, + V, X i3z (44) 

and 

^^P cR B B^V.x-jzxB,) • 

-^^ = ^^^i w^^ + 

The variational equation for x incorporates the vertical component of the induction 
equation (the part that depends on dB^/dt), but does not depend at all on d^h/dt. This 
means, for the variational solutions, that any observed time behavior for B/j can be specified 
independently of the time behavior for B^. 

To illustrate the variational technique, we minimize E"^ integrated over the magne- 
togram, {W'^ = 1) again using the ANMHD simulation described in §3.11 The "observed" 
map of dBz/dt is used as input, and we solved equation (H3|) for x, using Neumann boundary 
conditions (assuming zero horizontal Poynting fiux entering at the horizontal boundaries). 
This is a very good approximation for all but a few short segments of this synthetic magne- 
togram boundary. Since equation (H3i) shows that x = is a solution in regions of insignif- 
icant Bzi corresponding to the low field strength regions of the domain, the solution for x 
is set to zero for magnetic fields strengths below a threshold value. Once a solution for x is 
found, we verified that the horizontal components of E found from equation fl39|) obeyed the 
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Fig. 6. — Top panels show the values of the x— , y—, and z— components of cE taken 
from the ANMHD simulations described in the text. The middle panels show the same 
components of E determined with the variational formalism (equation [33], plus equations 
|D3] and [39], i.e. minimizing E"^). The bottom panels show scatter-plots of the original 
versus derived values of the x—, y—, and z— components of cE. All of the grayscale images 
are plotted using the same scale. 



- 24 - 



induction equation, i.e., c^h x E/j = —dBz/dt. The resulting three components of E are 
shown as the middle panels in Figure [61 

The results show a generally poor agreement with E from the ANMHD simulation. 
Thus, at least in this case, the variational approach of minimizing E"^ does not do a good 
job of reproducing the actual electric field. The PTD solution by itself {i.e. assuming that 
E = E^) shows better agreement with the simulation data, even though it produces spurious 
components of E parallel to B (which the ANMHD simulations did not have). Computing 
the rank-order correlation coefficients between the ANMHD and variational electric fields in 
the strong magnetic field regions results in values of 0.33, 0.50, and 0.50 for the x— , y— , 
and z— components of E, also indicating a poor agreement with the ANMHD electric field 
results. 

In spite of the poor comparison with the ANMHD results, the variational method did 
what it was designed to do: find an electric field that obeyed the induction equation, yet 
do this with minimum amplitude. The "halos" shown in the PTD solutions of Figure HJ for 
example, have been eliminated. Although the resulting electric field was not consistent with 
the simulation data, the variational method does yield a physically reasonable result, and 
stays zero in regions where one expects to find E near zero. These results motivate future 
exploration of other choices for W"^, to determine whether other choices result in better fits 
of the variational results to the MHD simulation data. 



4. An Example: NOAA AR 8210 



To illustrate the ideas described in the previous sections with an example using a real se- 
quence of vector magnetograms, we apply these techniques to a pair of vector magnetograms 
of NOAA Active Region 8210 taken with the University of Hawaii's I maging Vector Mag - 
netograph instrument (IVM) at Mees Solar Observatory on Haleakala (IMickev et al.lll996l ). 



These observations have already been described in detail in IWelsch et al 



f l2004h . We chose 



this observational test case because the data is well known, and we can therefore forego a 
detailed discussion of the data and its analysis. 

The pair of vector magnetograms, separated by a period of roughly 4 hours, were first 
cross-correlated and shifted to remove a mean shift due to solar rotation, and then averaged to 
define a mean vector magnetic field, and differenced to approximate a partial time derivative 
of each magnetic field component. 



We first apply the PTD formalism for the average vector magnetic field to derive the 
three fields i3, dB/dz, and J. In §2.31 we pointed out that the solutions are degenerate, in 
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that solutions to the homogeneous Cauchy-Riemann equations can be added to the solutions 
for OB/ dz and J without affecting the derived values of B/^. The degeneracy can be removed 
by choosing to set either dj jds = 0, or d/ds [dB/dz) = (but not both) when applying 
boundary conditions f[T7j) and f ITS]) to the solutions of the Poisson equations for these two 
functions. We show in Figure [2] how the functions dB/dz and differ depending on which 
parallel derivative is set to zero along the magnetogram boundary. Neither choice affects 
the derived values of B/^, but there is a slight advantage to choosing to set dJ /ds = 0: 
in that case, the observed values of i?„ are due entirely to d/dn{dB/dz), and therefore 
correspond to the potential-field solution that matches i?„ at the boundaries (Appendix [A|) . 
That means that any contribution to the horizontal magnetic field from currents can be 
attributed entirely to the contribution from J'. 

These points are illustrated in Figure U\ which shows the spatial distribution of B, 
Ap = X Bz, Bz, B/j, and the contributions of both the potential-field and current 
sources to B/j. To remove the Cauchy-Riemann degeneracy, it was assumed that dJ /ds = 
along the magnetogram boundary, coinciding with the choice of the top two panels of Figure 
[2J A homogeneous Neumann boundary condition is used to compute B: dB/dn = along 
the magnetogram boundaries. 

Considering next the time evolution of the magnetic field, we first apply the PTD 
solution to the measured values d^/dt. To solve equation for B, we have to assume a 
boundary condition for E at the edges of the magnetogram. Generally, this is not known, but 
if one is sufficiently lucky to have the magnetogram boundary located entirely in a region with 
zero or small magnetic fields, it is reasonable to assume that E = along the magnetogram 
boundary. If one is solving only the equation for B (and not solving for ip simultaneously), 
one can only require that one component of E/^ vanish along the boundary. We believe 
it is more physical to set the component of E/^ parallel to the boundary (Eg) to zero at 
the boundary, which means setting the normal derivative of B to zero there (homogenous 
Neumann boundary conditions), motivated by the discussion in §2.31 

Our data for AR8210 includes some regions that have significant magnetic field changes 
along short sections of the magnetogram boundary. Since we are merely trying to demon- 
strate how to use these techniques with real data, our approach here is to set dB/dt = 
in a narrow strip of three zones just inside the magnetogram boundary, and assume we can 
then set the parallel component of E/j (Es) in this slightly altered test case to zero. This 
boundary condition is equivalent to having a zero derivative of B in the direction normal 
to the boundary (Neumann boundary conditions). Once this assumption is made, it is also 
advisable to use the same Neumann boundary conditions for B, so that the time evolution 
of A and Ap are consistent with E computed at the boundaries. 
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If one is using only the PTD solutions, it is probably also a good idea to set a physically 
reasonable boundary condition for J at the edge of the magnetogram. For reasons similar 
to those described above, it seems reasonable to assume that = along magnetogram 
boundaries that lie in regions of weak average field. Since cE^ = — J in the PTD formalism, 
this can be achieved by choosing to set dj jds = when applying the boundary conditions 
(11711181) to solving the Poisson equations (jMHl) • Making this assumption removes the Cauchy- 
Riemann degeneracy for the JT" and dB/dz solutions. This has the added consequence of 
allowing one to interpret the contributions of dB/dz to d^h/dt as being changes to the 
potential-field part of the solution. 

Figure [8] shows the resulting PTD solutions for E^, and for V x E'^. Also shown are the 
decompositions of the curl into the evolution of the potential-field, and those driven by the 
observed evolution in J^. 

We have also tested the iteration technique and the variational technique on the AR 
8210 vector magnetogram data. Figure [9] shows a comparison of E computed using the three 
different techniques. The two left panels show E/^ and E^ from the PTD solution, with the 
boundary conditions as described above. The two middle panels show the same quantities 
from the iterative technique, and the two right panels show the solution from the variational 
method. 

There are common patterns in seen in all three solutions, namely that Eh swirls 
around a region of decreasing positive Bz on the right hand side of the magnetogram about 
40% of the way up from the bottom, and all three show a common pattern near the large 
sunspot. But the PTD solution also shows clear evidence of artifacts as well, such as a 
strong horizontal electric field normal to the magnetogram boundary near the left edge that 
is non-existent in the variational solution, and much less pronounced in the iteration solution. 
The iterative technique shows that some of the artifacts in the PTD solutions are reduced, 
but also shows large amplitude signals in some of the weak field regions. The variational 
solution shows small electric field vectors in regions of small magnetic field strength, which 
seems physical, at least superficially. There is little resemblance between the E^ solution 
found with PTD, and the E^ maps of the iterative and variational methods. Most likely, 
this is because the latter two methods enforce E ■ B = 0, meaning that E^ must be adjusted 
during solution algorithm such that this condition is satisfied. The iterative and variational 
solutions for E^ are poorly constrained near polarity inversion lines, judging by the large 
fluctuations near them, probably also a consequence of forcing E ■ B = 0. This might 
indicate that evolution is not consistent with ideal MHD along polarity inversion lines. 

Finally, we show the Poynting fluxes found from all three electric field techniques in 
Figure [TOl In contrast to the electric fields themselves, the Poynting flux distributions are 
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all fairly similar, showing a northward (positive in the y direction) flux of magnetic energy 
out of the negative polarity sunspot, and positive vertical Poynting flux in the region south 
of the negative sunspot. 



5. Summary and Discussion 

In this paper, we introduce and demonstrate three new techniques for determining elec- 
tric fleld distributions given a sequence of vector magnetogram observations. The following 
is a summary of the most important points. 

The flrst technique, based on a poloidal-toroidal decomposition (PTD) of the time 
derivative of B and hence cV x E, results in a solution for E that satisfles all three com- 
ponents of Faraday's law, but does not constrain any components of E derived from the 
gradient of a potential ip. We showed in §3.11 and in §l]that the unmodifled PTD solutions 
(assuming V?/' = 0) result in several artifacts and non-physical effects, and argued that more 
physically realistic solutions for E must include contributions to the electric fleld from a po- 
tential function, —'Vip. One concern in particular is that the PTD solution has signiflcant 
components of E parallel to B, contradicting the ideal MHD assumption that many regard 
as most likely to represent the solar photosphere, and certainly contradicting the ANMHD 
simulation test case discussed in detail throughout §21 In spite of these defects, the electric 
flelds found from the PTD solutions were more accurate in the strong magnetic fleld regions 
of the ANMHD sim ulation than those d etermined from all but one of the velocity inversion 



techniques tested in lWelsch et al.l (120071 ) . and were as accurate as those of the MEF method. 
The PTD solutions for E^ also enable easy-to-compute estimates of Poynting and helicity 
fluxes. 



We then explored two techniques for computing a potential function ip. First, in §3.2[ 
we outlined an iterative procedure for computing a solution —'Vip that acts to cancel the 
components of E parallel to B. This technique looks promising when applied to the ANMHD 
test case, signiflcantly reducing scatter of the inverted electric fleld values with those from 
the ANMHD simulations. When compared only in the strong fleld regions analyzed by 



Welsch et al.l ( 120071 ). the iteration method performed at a similar level to the PTD solutions. 
The improvement to the solutions seen in the scatter plots evidently occurs mainly in the 
weaker magnetic fleld regions. The solutions derived from this method are not unique, 
however, and when applied to the NOAA AR 8210 example (§!]), resulted in large electric 
fleld values in regions of small magnetic fleld, which would imply unphysically large velocities 
if an ideal MHD interpretation is assumed. 
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Fig. 7. — The upper left panel shows B as the background image, while Ap = V/^ x 
Bz, the vector potential for the potential magnetic field with the same Bz as the vector 
magnetogram, is shown as the arrows. The arrows in the upper right panel show B/i, which 
can be decomposed into the two contributions shown in the lower two panels. The lower 
left panel, also with as the background image, shows the contribution to B^ solely from 
V/i X J'z, which shows the contributions from non-zero values of J^. The lower right panel 
shows the vertical magnetic field Bz as the background image, and the horizontal components 
of the potential magnetic field = ^^dB/dz (equation |A4] ) as the arrows. All of the 
horizontal magnetic field vectors are drawn at the same scale. 
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Fig. 8. — The upper left panel shows cE^ = —j as the background image, while cE^ = 
— V/t X i3z, the horizontal electric field vector responsible for the time evolution of from 
the vector magnetogram difference, is shown as the arrows. The arrows in the upper right 
panel show d'Qh/dt, with dBz/dt as the background image. The lower left panel, also with 
dBz/dt as the background image, shows the contribution to d'Qh/dt solely from — V/^ x jTz, 
which shows the contributions from non-zero values of dJz/dt. The lower right panel shows 
dBz/dt as the background image, and the time evolution of the horizontal components of 
the potential magnetic field, dB^/dt = —'VhdB/dz as the arrows. The vectors in the three 
panels showing changes to B/j are all drawn at the same scale. 
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Fig. 9. — The top three panels show the derived horizontal electric field using the three 
different techniques discussed in this paper, applied to the AR 8210 vector magnetogram 
sequence. The background image in all three cases is dBz/dt. The arrows in the upper left 
panel show values of from the PTD algorithm with —'Vip set to zero. The top middle 
panel shows E/j = — ^hi', with ip computed using the iteration technique. The top right 
panel shows E/^ computed with the variational technique (minimizing E'^). The bottom three 
panels show E^^ (left), E!, — dip/dz computed with the iteration technique (middle), and Ez 
from the variational method (right). Note that Ez is poorly behaved in the latter two cases 
along magnetic neutral lines. All of the electric field vectors are drawn at the same scale, 
and the Ez images are displayed with the same linear grayscale colormap. 
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Fig. 10. — The top three panels show the derived horizontal Poynting flux using the different 
techniques discussed in this paper, applied to the AR 8210 vector magnetogram sequence. 
The background image in all three cases is Bz . The arrows in the upper left panel show values 
of from the PTD algorithm with —'Vip contribution set to zero. The top middle panel 
shows Sh computed from — ^ip using the iteration technique. The top right panel shows 
Sh computed with the variational technique (minimizing E"^). The bottom three panels show 
Sl (left), Sz computed with the iteration technique (middle), and 5*^ from the variational 
method (right). The vectors in the top three panels, and the grayscale images in the bottom 
three panels are all plotted with the same scale. For a better view of the Poynting flux 
structure, only a 60x60 pixel region of the vector magnetogram fleld of view is shown. 
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Second, in §3.31 we explored a variational approach to computing ip, by demanding that 
the total electric field minimize a positive definite integral over the magnetogram field of 
view, and that it obey the constraint E ■ B = (though the formalism allows for relaxation 
of this constraint) . This is essenti ally the same concept as Longcope's Minimum Energy Fit 



first, the minimization integral is generalized beyond the kinetic energy case considered by 
Longcope; second, we discovered that the variational solution for ip can be combined with 
E"^, resulting in a single equation for a new scalar function x whose gradient is proportional 
to the horizontal Poynting flux. The total electric field, including both the inductive and 
potential contributions, is then computed post facto from Vx- One can find V?/^, if desired, 
by subtracting the PTD solutions E^ from the full electric field solutions derived from Vx- 

We applied the variational technique to the ANMHD simulation results, and to the 
AR8210 vector magnetogram pair. The variational method does a poor job of reproducing 
the electric fields in the ANMHD simulation, at least in the one case we tried of minimizing 
E"^ {i.e., setting W"^ = 1). On the other hand, the solutions do correctly solve the induction 
equation, and they do so with a much smaller electric field than the PTD or iteration method 
solutions, or even the actual electric fields from the ANMHD simulation. The variational 
technique applied to AR8210 shows large electric fields mainly in regions where the magnetic 
field is changing rapidly, with the electric fields elsewhere being very small. Both the PTD 
and iterative solutions show large electric field vectors in regions of insignificant magnetic 
field. In this case, the variational solution seems more physical, though this is a subjective 
evaluation. 

We conclude that with only the magnetic field measurements, it simply is not possible to 
recover the true electric field in the solar photosphere - too much information is missing, since 
the true electric field depends not only on the induction equation, but also on solutions of 
the momentum and energy equations, about which vector magnetograms provide no direct 
information. But it is possible to find an electric field solution that is both physically 
reasonable, and consistent with the observed evolution of B. Our results thus far indicate 
that the variational technique does minimize electric field artifacts in regions where we don't 
expect significant electric fields. 

The PTD formalism for B leads to a useful and easy decomposition of the observed 
magnetic field into potential and non-potential contributions to the field. One appealing 
application of this result is that it suggests a recipe for evolving magnetograms as bottom 
boundaries of MHD simulations, from potential field distributions toward the actual magnetic 
field distribution. By starting from an initial potential field model in which = (see 
Appendix |A]), and then constructing a synthetic time evolution J' that evolves toward the 




However, we extend his technique in two important ways: 
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observationally determined distribution of J from the static PTD solutions, one could impose 
cE^ = — J (along with E/^ = 0) and allow an MHD simulation to evolve a solar atmospheric 
model toward the observed magnetogram state. It is important to emphasize, however, this 
is generally not consistent with an ideal MHD model for the photospheric electric field - 
one would need to have an MHD code with the flexibility to accommodate a user-specified 
electric field at the simulation boundary. 

The iteration and variational solutions, since they explicitly enforce E ■ B = 0, could 
be used both to define velocity fields at the photospheric bound ary, and for assirnilatin g 



a time series of vector magnetograms directly into MHD models (lAbbett fc Fisher! |2010| ). 
In the future, these solutions will be more thoroughly compared and contrasted with the 
more conventional inductive and tracking-based solutions for the velocity field from vector 
magnetogram data. 
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A. Potential Magnetic Fields Described with the PTD Formalism 

The PTD formalism allows for an alternative approach for describing potential magnetic 
fields near the magnetogram layer. Normally, potential magnetic fields derived from magne- 
tograms use the normal component of the field on the photospheric boundary, plus boundary 
conditions for the magnetic field at the side and upper boundaries (sometimes taken to be at 
infinity) to derive a solution within a specified volume. The horizontal fields of the potential 
solution on the bottom boundary are then determined by taking the horizontal gradient of 
the scalar potential that describes the potential field. The horizontal components of the 
potential field at the photosphere thus depend indirectly on the assumed behavior of the 
field at a significant distance from the photosphere. In this Appendix, we show how the 
horizontal field components in a vector magnetogram can be used to define a potential field 



solution using the PTD formalism, as an alternative to using an assumed behavior of the 
field at distant boundaries. 



Using equation ([T]) and setting all three components of the electric current to zero, one 
can show that V^JT" = and V| (V^B) = 0. This condition can be met with J' and V^B 
being functions of z only. We will assume henceforth that both functions of z are equal to 
the special case of zero: 

\/^B = J = (Al) 

Since a potential magnetic field can be specified with a single potential function, we assume 
we that we can find a single B function that can represent an arbitrary potential field, and 
ignore J'. 

What is the relationship between the usual scalar potential function (B = — V0) and 
the corresponding B function for the same potential field? To distinguish between the general 
case and the potential-field case, we denote the potential field case of B as Bp. Further, here 
we consider Bp to be an explicit function of three-dimensional space, in contrast to the 
convention in the rest of the paper that the scalar potentials are assumed to depend only 
on the two-dimensional domain of the vector magnetogram. From equation ([1]) without the 
contribution from J', 

B = V X V X Bpz . (A2) 

If the field is potential and thus current-free, then from equation ( 1A1[) Bp satisfies the three 
dimensional Laplace equation 

V^Bp = . (A3) 
In the PTD formalism, the horizontal field on the plane of the magnetogram is given by 

B. = V.^ (A4) 



while the vertical field is given by 



VlBp = -B, . (A5) 



However, since Bp obeys the three-dimensional Laplace equation, the left hand side of equa- 
tion (]A5|) is also equal to —d'^Bp/dz'^, and therefore, the magnetic field can be written as 



dz 

meaning that one can then make the identification 
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The PTD formalism allows one to use two different parts of the observed data in com- 
puting properties of the potential field. As noted above, the conventional approach uses the 
observed normal component [Bz) of the field at the magnetogram as the only photospheric 
boundary condition. The PTD formalism allows one to include the divergence of the hori- 
zontal field in the solution, as well as the normal component of B. Equation ([8]), applied to 
Bp, shows that 



This uses the horizontal magnetic field data to specify the rate at which Bz decreases im- 
mediately above the surface of the magnetogram in the potential field solution. To make 
the solution of the two-dimensional Poisson equation (lASIl well-posed, one can apply the 
Neumann boundary condition at the edge of the magnetogram. 



where as in §2.11 Bn is the observed component of B/j normal to the magnetogram boundary. 

From the solution to the horizontal Poisson equation (IA8|) for dBp/dz, one can then 
use equation (IA4p to reconstruct the contribution to B/i that can be ascribed solely to a 
potential magnetic field, without having to assume any behavior at distant boundaries. On 
the other hand, using the potential field derived in this way to compute the behavior far 
from the magnetogram is probably dangerous, as any errors in the field measurement are 
likely to be greatly magnified when extrapolated to large distances. 

Finally, it is frequently useful to be able to express the potential field in terms of a vector 
potential Ap, rather than a scalar potential. Estimates of the magnetic helicity flux through 
the photosphere typically involve Ap. Using the PTD formalism, the vector potential Ap 
for the potential field is found simply from equation ([2]): 



evaluated in the plane of the magnetogram, and where Bp is found from a solution to the 
two-dimensional Poisson equation (IA5p . Note that Ap has no vertical {z) components. 




(A8) 




(A9) 



Ap = V/i X Bpz , 



(AlO) 



B. Fourier Transform Solutions to the PTD equations 



If the boundary conditions for the vector magnetogram are assumed to be periodic, and 
if the magnetic field's time derivative does not have a uniform (zero-wavenumber) compo- 
nent. Fast Fourier Transform (FFT) techniques greatly simplify the solutions for the Poisson 
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equations for B, J7, and dB/dz. If we denote the Fourier transforms of B^, By, and Bz as 
Bx, By, and Bz, respectively, one can write the solutions to equations (EHH]) as 



and 



(B3) 



where J^~^ denotes the inverse Fourier transform, and and ky are the horizontal wavenum- 
bers in the x— and y— directions within the magnetogram, respectively. The quantities B 
and J can then be used to derive the inductive electric field via equation ^\M . 

The same approach can be used to find the functions B, J , and dB/dz from Fourier 
transforms of the magnetic field B^, By, and Bz by using the same equations flBmB3P except 
with the Fourier transforms of the time derivatives of the magnetic field components replaced 
with the Fourier transforms of the magnetic field components themselves. 

However, if the magnetic field, or its time derivative, has a non-zero average value in 
any of the component directions, the FFT solutions cannot account for this, and one must 
use the techniques described in Appendix [C] to correct the FFT-derived solution. Equation 
02]) can then be used to to find the vector potential (ignoring the gauge contribution). 



C. Accounting For Average Values of the Magnetic Field and its Time 

Derivative 

If the magnetic field time derivative contains a spatially uniform (zero-wavenumber) 
component Bq, equations (lBltlB2p will not recover this component, because the assumption 
of periodic boundary conditions for E cannot produce a uniform vector V x E and hence B. 
An additional electric field component must be explicitly added to account for the uniform 
time derivative of B. Similarly, if Neumann boundary conditions are used to solve equation 
dnj, the solution will force an assumption that the spatial average of Bz is zero, and an 
additional term must be added to the derived solution for E. 



A term of the form ^ 

cEo = --Bo X (r - ro). 



(CI) 
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will fully reproduce the observed time derivative of the magnetic field when added to the 
electric field derived by using equations f lBltlB2p in equation (JHj). Here r is the position 
vector, and tq can be any constant vector offset. 

To correct equation ( IT3|) for cV x E, it is sufficient to just add the term — Bq to equation 
(fT3|) if using the FFT-derived solutions for J and dB/dz given in Appendix [Bl 

A similar problem arises if the magnetic field itself has a non-zero uniform component 
Bq: the Fourier transform-derived solutions for A will not be able to recover Bq. Instead, a 
term of the form ^ 

Ao = - Bo X (r - ro) (C2) 

when added to equation ([2]) will reproduce the full magnetic field observation as V x A. 
Note that Eq = -1/c dAo/dt. 

The solutions for cEq depend on an unspecified position vector offset, tq. We now 
present an argument for determining tq, based on the concept of Galilean invariance. 

If the electric field contribution from equation (IC1|1 originates from an ideal electric field 
due to a velocity field vq in the presence of the uniform magnetic field component Bq, then 
we can equate the two expressions for the electric field: 

-lBox(r-ro) = -^xBo. (C3) 

Taking the cross-product of both sides of equation flC3l) with the vector Bq, one finds after 
some manipulation 

i52(r-ro)^ = (voxBo)xBo, (C4) 

where Bq is the amplitude of Bq, and subscript _L denotes the directions perpendicular to 
Bq. Taking a spatial average of equation flC4p then results in 

^52(r - ro)^ = (Vo X Bq) x Bq. (C5) 

Now we can identify the spatial average vq with Vj-e/, the velocity of a Galilean reference 
frame. In other words, if two magnetograms in a sequence have an overall non-zero net shift, 
due to e.g. an inaccurate (or non-existent) correction for solar rotation, we then identify the 
overall reference frame velocity responsible for this shift as v^e/. If the two magnetograms 
have been co-registered such that the net overall shift is zero, then we assume that v^e/ = 0. 
In that case, the right hand side of equation ( \C5\i must be zero, and one then finds this 
constraint for tq: 

ro± = f^. (C6) 
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This condition is always satisfied when the vector offset fq, which can also be regarded 
as the origin of the vector magnetogram coordinate system, is chosen to coincide with the 
geometric center of the vector magnetogram f . By determining Tq, this allows us to determine 
an electric field and vector potential solution unambiguously. If one wants to perform the 
calculation in a moving reference frame with a non-zero v^gj, the value of fq will then be 
determined from equation ( 1C5|) instead of equation ( ]C6p . 

It is useful to evaluate the result after applying this substitution into equations (ICip 
and (lC2p in component form for a vector magnetogram lying on the surface z = 0. In that 
case, X and y represent the average values of x and y at the center of the magnetogram, and 
^ = 0. We find 

cEo = ^[{y- y)Bo^z -{z- z)BoJx 

+i[(x-x)4,o- (z/-y)4,o]z , (C7) 

and 

Ao = ^[-(y - y)Bo,z + z)Bo.y]^ 
- x)Bo^^ -{z- z)Bo^^]y 
+ ^[-{x-x)By^o + {y-y)B,,o]^ , (C8) 

where -Bq.x, Bo^y, and Sq.z are the components of Bq, and -Bo,x, Bo^y, and -Bq.z are the 
components of Bq. While the terms proportional to {z — z) will vanish at the photosphere 
in Aq and cEq, they will still contribute to vertical derivatives of Aq and cEq evaluated at 
the photosphere. 



D. Deriving the Variational Equation for x 

Equation ( 137|) can be re-written as 

{W^/Bz){cY.h X BzT. + cEzT. X B^) = -VnX. (Dl) 

or 

cE;, X z + {cEz/Bz)z xBh = — (D2) 
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Taking the dot product of equation (ID2I) with z x B^, and making use of equation ([3 
namely -Ez-B^ = (R ■ B — E/j ■ B/j), we find 

= Bf-^ - . z X B,.. (D3) 

Substituting equation flD3P into equation (ID2p then results, after some manipulation, in 

cE, X z + z X B, + :^(B. ■ V,x)B. = "^^V.x (D4) 

At this point, the only quantity not involving x that is unknown in equation ( ]D4p is 'Eh 
(recall that R is assumed to be either zero or known a priori). If we take the horizontal 
divergence of equation ( ]D4I) . however, we can eliminate cE/j through the use of the magnetic 
induction equation V/i ■ (cE/^ x z) = cz ■ V/^ x E/i = —dBz/dt resulting in 

Qj^ /cR B \ / 1 \ / \ 

-df + ■ \^ z X B,j + V, ■ [y^.iBh . V,x)B.j = -V, ■ [y^.^kXj ■ 

(D5) 
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